Topological phases in a two-dimensional lattice: Magnetic field versus spin-orbit coupling 
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In this work, we explore the rich variety of two-dimensional topological phases that arise when considering 
the competing effects of spin-orbit couplings, Zeeman splitting and uniform magnetic fields. We investigate 
minimal models, defined on a honeycomb lattice, which clarify the topological phases stemming from the in- 
trinsic and Rashba spin-orbit couplings, and also from the Zeeman splitting. In this sense, our work provides an 
interesting path connecting the quantum Hall phases, generally produced by the uniform magnetic field, and the 
quantum spin Hall phases resulting from spin-dependent couplings. First, we analyze the properties of each cou- 
pling term individually and we point out their similarities and differences. Secondly, we investigate the subtle 
competitions that arise when these effects are combined. We finally explore the various possible experimental 
realizations of our model. 
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I. INTRODUCTION 

Today, unprecedented efforts are devoted to the study of 
novel topological insulating phases, which exhibit remark- 
able properties at their boundaries.' ^ In two dimensions, these 
quantum states are characterized by propagating edge states, 
which carry dissipationless currents along one-dimensional 
boundaries. From a theoretical point of view, these trans- 
port properties rely on the existence of large bulk gaps in the 
energy spectrum, which host robust and gapless edge excita- 
tions. Different edge state structures and transport properties 
can be realized according to the nature of the gaps: A mag- 
netic field breaks time-reversal symmetry (TRS) and leads to 
chiral edge states (i.e., paiticles propagating in a given di- 
rection) and quantized Hall currents,^ whereas a large spin- 
orbit coupling preserves the TRS and produces helical edge 
structures. ^ In the latter situation, referred to as the quantum 
spin Hall (QSH) effect, edge states with opposite spins coun- 
terpropagate and contribute to a vanishing Hall current, while 
producing a spin Hall current.^ ^ The origins of the quantum 
Hall (QH) and QSH phases are fundamentally different, as 
the former is produced by an external field, whereas the lat- 
ter relies on the intrinsic properties of the material. However, 
one can interpret the QSH phase as being two opposite QH 
phases, one for each spin component. This observation is eas- 
ily shown through the Kane-Mele lattice model,^ ^ which real- 
izes the QSH effect by the inclusion of the intrinsic spin-orbit 
coupling (ISO), and corresponds to two copies of the (QH) 
Haldane modeL- 

In recent works, the ISO coupling has been combined with 
a magnetic field,^ '" or with a constant exchange term,^in or- 
der to study the effect of the TRS breaking on the QSH effect. 
It has been shown that the QSH effect persists in the presence 
of magnetic field, leading to the "TRS -broken'" or "weak" 
QSH effect.- The term "weak" refers to the absence of robust- 
ness against spin-flip scattering due to magnetic disorder. In 
absence of this scattering process, the Hall and spin Hall con- 
ductivities are protected and quantized. 

In this work, we present a different approach to study the 



interplay between QH and QSH physics. Contrarily to the 
general trend in the field of topological insulators, which is to 
include many terms and study complicated Hamiltonians, here 
we investigate the minimal models that can produce topolog- 
ical phases. We consider a tight-binding model of the hon- 
eycomb lattice under a perpendicular magnetic field, which 
includes the ISO coupling, the Rashba spin-orbit (RSO) cou- 
pling, and the Zeeman effect. This spin- 1/2 model was intro- 
duced in our recent paper,l21 where we already studied the com- 
bined effects of some of these terms. Here, instead, we first 
investigate the effect of each of these terms individually on the 
QH phases generated by the magnetic field, in order to distill 
the problem. In particular, we demonstrate the equivalence of 
the ISO and Zeeman terms in generating the weak QSH phase 
at zero energy in the presence of a magnetic field. Then, we 
show that the ISO coupling can drive topological phase transi- 
tions between different topological insulating phases, an effect 
which totally relies on the combination of the ISO coupling 
and the magnetic field. 

Secondly, we study the interplay between the terms incor- 
porated in our model. For example, we discuss the competi- 
tion between the ISO coupling and Zeeman effect, motivated 
by the similarity of these terms for generating the weak QSH 
effect. We show that the combination of these terms does not 
necessarily lead to a larger regime of parameters where the 
weak QSH effect is observed, which reveals the subtle com- 
petition between each of these terms and the magnetic field. 
We also study the exquisite effects produced by the Rashba 
coupling in the presence of a large exchange field, which is 
shown to generate several QH gaps, even in the absence of a 
magnetic field. From another perspective, we also discuss the 
fate of the weak QSH phase, as the RSO coupling is turned on, 
which introduces spin-flip terms into the Hamiltonian, so that 
the spins are no longer aligned perpendicularly to the plane, 
and the spin Hall conductivity is no longer protected. 

An important motivation for the investigation of these topo- 
logical phases and phase transitions is the rich variety of sys- 
tems where they may be observed. Despite its weak ISO cou- 
pling, graphene remains an appealing candidate for the obser- 
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vation of these phases, e.g., by increasing the ISO coupling ar- 
tificially by adatoms.'^ Secondly, several experimental groups 
have reported on synthetic honeycomb lattices in a condensed 
matter system, for instance arrays of quantum dots on a GaAs 
substrate,'^ and "molecular graphene", where the lattice is 
created by using the repulsive interactions of CO molecules 
deposited in a triangular array on a Cu(lll) surface."^ These 
systems have larger lattice constants, which lead to a higher 
flux per plaquette at realistic magnetic field values, and al- 
low for a more flexible control over the parameters than in 
real graphene. A third type of experiments that motivates our 
studies is the recent realization of a honeycomb optical lattice 
for ultracold atoms. '^ '^ In cold-atom systems, large magnetic 
fieldJ^and spin-orbit couplings''' are produced synthetically 
by adjustable lasers (see Ref. 20 for a review). In particu- 
lar, such artificial gauge potentials could be produced in op- 
tical lattices,'22lEIl thus leading to the possibility to probe the 
high-flux regime. Furthermore, we note that optical square 
lattices subjected to well-designed gauge potentials can repro- 
duce the properties of honeycomb lattice systems, exhibiting 
Dirac-type physics.'^—' Cold-atoms systems carry the advan- 
tages of high tunability, necessary for studying phase transi- 
tions, and allow for the study of topological phases in absence 
of interparticle interactions^*^ and disorder.I^Sl These possible 
experimental realizations are explored in detail in the final part 
of this work. 

This article is structured as follows. We introduce the lat- 
tice model in Sec.|ll]and present the four terms characterizing 
our tight-binding Hamiltonian: the usual hopping term, the 
two spin-orbit terms as well as the Zeeman splitting. Then, 



in Sec. Ill we review the techniques used for our investiga- 
tion and characterization of topological phases. Section [rV| 
discusses the combined effect of the magnetic field and spin- 
orbit coupling. Then, the combined effects of spin-orbit cou- 
plings, Zeeman splitting and uniform magnetic flux are ex- 
plored in Sec.[V] The possible experimental realizations of our 
model, as well as the possibility to detect the effects presented 
in Sees. IV and|V] are reported in Sec. VI Finally, the conclu- 
sions are drawn in Sec. VII The Appendix provides further 
mathematical details of the calculations. 



II. MODEL 

We consider a tight-binding model of spinful electrons in a 
two-dimensional honeycomb lattice, subjected to a uniform 
perpendicular magnetic field B ~ Be^- This model was 
introduced in Ref. |9] to investigate the effects of an external 
magnetic field on the Kane-Mele model (cf. Refs. 4 and 5). 
The magnetic field induces a Zeeman effect, and causes the 
hopping terms to acquire a phase. This phase is encoded by 
the so-called Peierls substitution, where one replaces the mo- 
mentum p by p — eA, where A is the gauge potential as- 
sociated with the magnetic field. Thus, any hopping term 
from site k to site j picks up the phase factor e'^^*", where 



(e/H) A • dl. In the remainder of this text, we will 



of the lattice, expressed in units of the elementary flux quan- 
tum h/e. 

In our model, we include four effects into the Hamiltonian, 
written as H = Hnn + Hz + Hi + Hn, with 



-f^NN 

Hz 
Hi 

Hk 
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Here, Ck = (cfc^,Cfe^) is the annihilation operator for elec- 
trons at site k. The first term iJNN describes hopping between 
two nearest-neighbor (NN) sites j and fc, where t denotes the 
hopping amplitude. The second term Hz describes the Zee- 
man effect, which is an on-site term that assigns different po- 
tentials to the two spin components via the Pauli matrix cr^. Its 
amplitude is proportional to the magnetic flux (p and the co- 
efficient Az, which is related to the material's Lande g factor 
by 27r0Az = .<7Mb-B, where /ie is the Bohr magneton. Note 
that in a cold-atom emulation of this model, the parameters 



express the magnetic field strength (flux density) in terms of 
the dimensionless quantity 0, defined as the flux per unit cell 



(f> and Az are tuned individually (cf. Sec. VIC i. The Zeeman 
term lifts the degeneracy by shifting the two spin-degenerate 
copies of the spectrum up and down by an equal amount of 
energy. The third term Hi describes the ISO coupling, which 
corresponds to a next-nearest-neighbor (NNN) hopping with 
amplitude ti.* Here, the factor Vjk ~ ±1, where the sign de- 
pends on the value of d^j A dij, i.e., the outer product of the 
two bond vectors connecting site k to site j via their unique 
common neighbor at site /. This hopping also involves the 
matrix cr^, so that the sign of the hopping amplitude is oppo- 
site for the two spin components '[ and |. Thus, the two dif- 
ferent spin species are effectively subjected to opposite local 
(Haldane-type^) magnetic flux. In absence of a magnetic field, 
this term opens a topologically nontrivial gap, and causes the 
QSH effect, characterized by so-called helical edge states, 
which are protected by TRS.'*^ The final term H^ is the con- 
tribution to the NN hopping due to the RSO coupling. This 
hopping has an amplitude tn, and involves the spin matrix 
(7 A djk = cFxd^jk — '^yd^k^ where cr — {a^, cry), and djk is 
the vector connecting sites k and j. The resulting spin matrix 
has only off-diagonal elements, so that this hopping involves 
a spin flip. This term couples the two spin components, and 
as a consequence, is no longer a good quantum number 
for tlx, 7^ 0. In practice, this means that spin states are no 
longer exclusively up or down, but may also point in differ- 
ent directions. Finally, we remark that the model considered 
here preserves inversion symmetry. In particular, we disre- 
gard the effect of an additional staggered potential, which was 
shown to induce topological phase transitions between trivial 
and non-trivial phases in Refs.mandjSj The possible physical 
realizations of our model ([T|) are discussed later in Sec.|VI| 
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III. TOOLS FOR THE TOPOLOGICAL ANALYSIS 

In the following sections, we will perform an analysis of the 
various topological phases produced by the terms in the model 
Hamiltonian Q. Before going into details, describing each 
term individually and their combined effects, we first present 
the framework and tools which we will use for this analysis. 
Readers who are already familiar with the notions of bulk and 
edge states, Chern numbers, and (spin) Hall conductivity may 
choose to skip this section. 



A. Harper equation and the Hofstadter butterfly 

We will illustrate our framework by studying the spin- 
degenerate model, i.e., we set Xz ~ h — tji = in Hamil- 
tonian ([T]). We assume t > throughout the whole text. The 
bulk band structure can be computed numerically by apply- 
ing periodic boundary conditions, namely by considering a 
toroidal geometry. This requires to set (p — p/q, where p and 
q are integers, in which case the system is periodic in both 
spatial directions. Under these conditions, and by choosing a 
proper gauge for the Peierls phases e*^J'=, the system reduces 
to a q X 1 magnetic unit cell. (See the Appendix for a descrip- 
tion of the gauge structure and its spatial periodicity.) By ap- 
plying a Fourier transform and invoking Bloch's theorem, the 
Schrodinger equation reduces to a 4(7 x eigenvalue problem, 
known as the discrete Harper equation,^" ^' 



7^L 1*. 



(2) 



where = (V'n At; V'n a;, V'nBt, V'nB;) dcnotcs the 

single-particle wave function at site n = 1, . . . ,q, and where 
the 4x4 matrices I?„ and Tin are given in terms of the pa- 
rameters defined in Hamiltonian ([T]l (see the Appendix for de- 
tails). The four components of the wave function ^„ are due 
to the two sublattice iA,B) and the two spin (t,i) degrees of 
freedom. The Aq x 4q Harper problem in Eq. Q is analogous 
to the original (spinless, square-lattice) Hofstadter problem,!^ 
where the Harper equation involves aq x q matrix. 

Solving the Harper equation (|2]i provides the bulk energy 
bands E = E{k.j;,ky). It is instructive to evaluate the 
sizes of the bulk gaps as a function of the magnetic flux <f>, 
which r esults in a diagram known as the Rammal-Hofstadter 
butterfly,'22l22l shown in Fig. [l] This figure may be obtained in 
the limit of large q, where the bands become dispersionless. 
The fractal structure, which is periodic in (j) with a period of 
1, is a result of the competition between two length scales: 
The magnetic length Ib = y/h/eB and the lattice spacing. 
Here, (p parametrizes this competition, as the ratio between 



the area of the unit cell and that of the cyclotron orbits ('■ 



B. Topological invariants 



1%)- 



Each bulk gap in the spectrum shown in Fig.[T]may be char- 
acterized by a topological invariant, which encodes the Hall 
conductivity of this system when the Fermi energy lies inside 
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FIG. 1. (Color online) Hofstadter butterfly spectrum for the spin- 
degenerate model, Hamiltonian |TJ with Xz ~ ti = in ~ 0. The 
bulk bands, i.e., the regions where the density of bulk states is pos- 
itive, are shown in black. In the bulk gaps, the colors and numbers 
indicate the quantized Hall conductivity in units of /h. 




FIG. 2. (Color online) (a) The cylinder geometry representing a rib- 
bon with periodic boundary conditions in the x direction, (b) Spec- 
trum at(j> — 1/3 for the spin degenerate model (Az = ti = — 0). 
The colors represent the expectation value (y): Red and blue indi- 
cate edge states at the bottom and top edge, respectively, and gray 
represents bulk states. The spin-degenerate edge states shown in (a) 
correspond to those in the QH gap at the Fermi energy indicated by 
the dashed line in (b). 



this gap. The integer values of the Hall conductivity are robust 
against external perturbations: They remain constant as long 
as the bulk gaps remain open. In the following, we recall the 
relation between this transport coefficient and the concept of 
edge states and topological invariants. 



1. Edge-state analysis 

The edge-state analysis can be performed for a system fea- 
turing a boundary,^^ such as a cylinder, as shown in Fig.|2|a). 
When solving the Harper equation (|2| on a cylinder, namely 
by applying periodic boundary conditions along one spatial 
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direction only, the spectrum consists of bulk bands and topo- 
logical edge states?^ Indeed, we typically find a few edge 
states within the bulk gaps, some of which cross the gap from 
one bulk band to the other; see Fig. |2|b). Importantly, each 
edge state contributes e'^/h (one quantum of conductivity) to 
the Hall conductivity of the system. The total Hall conductiv- 
ity is therefore equal to an integer number (i.e., the number of 
edge-state branches) times e^/h. The topological invariance 
of the Hall conductivity is due to the fact that the number of 
edge states inside the bulk gap cannot be altered unless the gap 
is closed. This result can be used to derive the topological in- 
variant associated to any of the bulk gaps, leading to the colors 
in the spectrum in Fig.[T] For this spin-degenerate model, all 
topological invariants are multiples of 2 because the system is 
twofold spin degenerate. 

More precisely, in order to evaluate the Hall conductivity of 
a bulk gap, we count the edge states whose dispersions inter- 
sect the Fermi energy, taking into account their location and 
direction of propagation. The location of each state is derived 
by computing its expectation value (y) from the eigenstate, 
which provides the coloring of Fig. |2|b). Secondly, the di- 
rection of propagation can be obtained from the sign of its 
momentum derivative dE/dk, where k is the momentum par- 
allel to the edge. Edge states with opposite directions con- 
tribute with opposite signs to the Hall current and therefore to 
the Hall conductivity, 

(TH = (iVR-iVL)^, (3) 

where A'l and Nj^ denote the number of left- and right- 
moving states, respectively. The Hall conductivity does not 
depend on the magnitude of the velocity (cx dE/dk), but in- 
stead is given by the number of edge-state channels. In this 
spin-degenerate model, the states are all either left- or right- 
moving, i.e., they are chiral states. The chirality of the edge 
modes provides protection of the Hall conductivity against 
backscattering due to disorder in the presence of broken TRS. 

2. Bulk topological invariants 



This identity is known as the bulk-boundary 
correspondence,^^ as it relates the topological indices 
associated to the bulk bands to the number of edge excita- 
tions. This correspondence provides us with an important 
observation: Although the dispersions of the edge states 
depend on the shape of the edges (e.g., zigzag, bearded, or 
armchair) and on the system size, the number of edge-state 
branches does not: this number is a bulk property and 
therefore does not depend on the specific form of the edge. 
In the remainder of this paper, we will study systems with 
zigzag and bearded edges, without any loss of generality 
concerning the topological properties. 

In the next sections, we will characterize the topological 
phases by their associated Hall conductivity and edge-state 
structures. As described in the following, the edge states pro- 
duced by the magnetic field are affected in the presence of 
strong perturbing terms, such as the spin orbit couplings and 
Zeeman splitting. For the sake of simplicity, we will show 
these effects using the edge-state analysis obtained from a 
cylindrical geometry, rather than focusing on the analysis of 
the bulk topological invariants. 



3. The Stfeda formula 

An efficient method to compute the Hall conductivity can 
be formulated in terms of the integrated density of states. In- 
side a bulk gap, the integrated density of states {E) is de- 
fined as the fraction of all states that lie below the gap (which 
is an integer times 1/Aq). The Hall conductivity is then pro- 
portional to the (/)-derivative of this quantity. 



a result known as the Stfeda formula.EHEni 

We have used this 

powerful formula to compute the Hall conductivity, and gen- 
erate the colorful butterfly spectra presented in the remainder 
of this paper. 



Remarkably, the number of edge states is directly related 
to topological invariants, the Chern numbers, a quantity asso- 
ciated with each of the bulk bands. The Chern number is an 
integer topological index defined in the toroidal geometryj ^^l^^ l 



1 



BZ 



dkx 



dk,, 



(4) 



where A„.k = (ti„.k| Vk|wri.k) is the Berry connection as- 
sociated to the eigenstate |u„,k), and where BZ denotes the 
(toroidal) Brillouin zone. The number of edge states Nr in- 
side the rth bulk gap, or equivalently the Hall conductivity 
assigned to this gap (cf discussion above), is then equal to 
the sum of all the Chern numbers associated to the occupied 
bands,tiE2l 



E 

occupied bands 



(5) 



4. Charge and spin Hall conductivities 

In the presence of terms involving the matrix cr^ in the 
Hamiltonian, the spin SU(2) symmetry is broken. If spin-flip 
terms (i.e., linear combinations of and ay) are absent, the 
spin in the z direction is conserved, i.e., it is a good quan- 
tum number. In this case, which happens if the ISO cou- 
pling or the Zeeman splitting is present while the RSO is not, 
the Hamiltonian may be decomposed into the two spin com- 
ponents (t,4-)- Thus, one is able to define the component- 
wise Chern numbers „ and Cj, „ and Hall conductivities 



(THt aiid ch^. The charge Hall conductivity of the whole 
system is the sum of the two component-wise conductivities, 
cth = cth ^ + cth ^, since both of them contribute equally to the 
Hall current. Considering their contributions to the spin cur- 
rent, spin-up and spin-down edge modes have opposite "spin 
charges" ±fi/2. Thus, the spin Hall conductivity is equal 
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to the difference of the two component-wise conductivities, 
cr^ /{e/Air) = (crnt ~ o'H^)/(e^//i), where e/A-K is the ele- 
mentary quantum of spin Hall conductivity. 

In the spin-degenerate model discussed earlier, the two 
component-wise Hall conductivities are always equal, so that 
all spin Hall conductivities vanish. This is certainly different if 
the spin degeneracy is broken. For instance, in the presence of 
ISO, the QSH gap at zero magnetic field and at zero energy^is 
characterized by ant = ^'^Hi = —e^/h. Here, the contribu- 
tions to the charge Hall conductivity cancel, but the contribu- 
tions to the spin Hall conductivity add up to cr^ = — 2e/47r. 

If the Hamiltonian contains spin-flip terms, e.g., the Rashba 
coupling, then the spin is not conserved, and the aforemen- 
tioned definition of the spin Hall conductivity is no longer 
valid. However, the related spin Chern numbers, which 
equal = „ — „ if spin is conserved, remain 

well-defined topo logical invariants, even in absence of spin 
cons ervation .EEl 



IV. EFFECTS OF THE MAGNETIC FIELD 

A. The spin-degenerate model: Landau Levels and the 
anomalous quantum Hall effect 

Zooming in on the butterfly of Fig. [T] near w 0, we ob- 
serve that the thickness of the bands decreases in the low- 
flux regime (not shown). Thus, in this limit we can treat 
the bands as being infinitely thin, i.e., as Landau levels. In 
the low-energy regime, between the Van Hove singularities at 
E/t = ±1, the band structure at low flux coincides with the 
Landau-level spectrum known from graphene,"^^ with ener- 
gies E cx ±\/2Til(f) (where Z = 0, 1, 2, ... is the Landau-level 
index). 

Moreover, computing the Hall conductivity in this region, 
we find the "anomalous" Hall-conductivity sequence an = 
4(n + l/2)e^/ft, (n E Z), as observed in graphene."*^ The 
breakdown of this structure around E/t = ±1 can be under- 
stood from the large density of states near the Van Hove singu- 
larities and from the Chern numbers associated to these many 
bands."**' The edge-state structures. Hall conductivity plateaus 
and topological aspects of spinless electrons in a honeycomb 
lattice subjected to a magnetic field have been thoroughly de- 
scribed in Ref..46i 



B. Intrinsic spin-orbit coupling 

Since the seminal work of Kane and Mele,''it is known that 
the ISO coupling opens a topologically nontrivial gap of size 
A — 6\/3ii at zero magnetic field. This bulk gap hosts two 
counterpropagating edge modes per edge, with opposite spins 
(one Kramers pair): These helical edge states are related by 
TRS. This topological phase is known as the QSH state, and 
may be regarded as two opposite QH phases (i.e., each spin 
performs the QH effect, with opposite chirality). The system 
preserves TRS (due to absence of the magnetic field), which 
protects the QSH state against scattering processes caused by 
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FIG. 3. (Color online) Hofstadter butterfly spectrum for the system 
with ISO coupling (ti/t = 0.1). Although the periodicity in is 6 
rather than 1 (see the Appendix), we present only the range < (p < 
1 for the sake of comparison with the other terms. The colors and 
numbers indicate the quantized Hall conductivity in units of /h. 



disorder. Although for any edge state a reversely propagat- 
ing mode would be available to scatter to, the spin has to be 
flipped, and the two different ways to flip the spin interfere 
destructively with each other, due to a mutual phase factor of 

The combination of the ISO coupling and a magnetic field 
leads to the breaking of the TRS satisfied by the QSH phase. 
However, due to the absence of spin flip terms in the Hamil- 
tonian, the helical edge-state structure persists. In terms of 
charge and spin Hall conductivity, the resulting state is equiv- 
alent to the QSH state. As announced in the introduction, 
this state is called the weak QSH phase^ or TRS broken QSH 
phase.E^From the Hofstadter butterfly in Fig. [s] we observe 
that the size of the weak QSH gap at zero energy tends to 
decrease if the flux is increased (but not monotonically), and 
closes eventually. In Fig.[3j in which tj = O.lt, the QSH gap 
at £^ = remains open in the range < < 0.2. Here we ob- 
serve the competition between the ISO coupling, which opens 
the (weak) QSH gap, and the magnetic field, which "tries" 
to destroy the weak QSH state by its TRS breaking property. 
Throughout the weak QSH gap, the edge states remain cross- 
ing at zero energy, which indicates that this phase is robust 
at least in the absence of magnetic disorder We finally note 
that the spectrum remains particle-hole symmetric, as in the 
spinless case. 



C. Zeeman effect 

The Zeeman term models an on-site spin-splitting effect 
that is always present in real materials subjected to mag- 
netic fields. Studying Hamiltonian ([T]i with Az 7^ and 
t]x = ti = 0, we observe that the spin degeneracy is lifted 
by merely shifting the two spin states up and down in energy 
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FIG. 4. (Color online) Hofstadter butterfly spectrum for the model 
described by Hamiltonian ([TJ with Zeeman effect, Xz/t = 0.5 and 
ti = tR = 0. The bulk bands are in black. The colors and 
numbers indicate the quantized Hall conductivity in units of /h. 
The diagonal dashed lines indicate the outline of the butterfly (i.e., 
E/t = ±3 ± 2-n4>\z/t) for the two spin components. 



by 2it\z4>- The spectrum is otherwise left invariant. With 
respect to the Hofstadter butterfly represented in Fig. [T] this 
means that the spin-up and spin-down copies are "skewed" 
in opposite directions, as shown in Fig. |4] Because the low- 
est and highest bands of the spin-degenerate model are situ- 
ated at approximately E/t = ±3, the two spin components 
are completely separated above a certain flux value, approxi- 
mately equal to = 3/(27rAz/t). Furthermore, the Zeeman 
shift breaks the periodicity of the butterfly spectrum along the 
(j) axis. 



up and down, as shown in Fig.|5|b) for = 1/61. Assuming 
that the splitting is sufficiently large, a new gap opens at zero 
energy, where the original band was situated. The spin-up and 
spin-down edge states now connect to bulk bands that have 
been shifted up and down, respectively. As a consequence, in 
the newly opened gap at zero energy, we have A^i^ = — 1 and 
— +1. Thus, the charge Hall conductivity vanishes, and 
the spin Hall conductivity is equal to cr^ — — 2e/47r. Based 
on these values, we find that this state is a weak QSH phase. 

Comparing this state to the one generated by ISO (see 
Fig. |5jc)), we observe the similarity between them. In both 
cases, the Landau levels for the two spin components shift in 
opposite directions. The difference between the two terms lies 
in the different values of the shift: The Zeeman shift ±27r(/)Az 
is linear in the flux, so that it vanishes for ((> ^ {), while 
the Landau-level shift induced by the ISO coupling equals 
zbSA/S^i.'^ A recent study shows that a flux-independent Zee- 
man term (called an exchange term) may also generate the 
QSH state at zero magnetic field. ' ' 

The similarity between the Zeeman and ISO terms may be 
understood at a formal level, in terms of the respective lin- 
earized (low-energy) Hamiltonians: The linearized Zeeman 
term acts as *'f((T|^ (g) 11^^ (g) ]1^^')\1/, where the factors 
indicate the proper spin, sublattice pseudospin, and valley 
pseudospin, respectively. On the other hand, the ISO cou- 
pling acts as \[''''(cr|^ ® af^ ® a^^')^. A priori, these 
Hamiltonians act differently on the eight-component field vec- 
tor \l/ = {iprT,T,i)a=^i,T=AB,£,=KK' ■ In the lowest Landau 
level, the four components V't.A.AT, ipi,A,K,^t,B,K' , ^i,B,K' 
vanish.'*'* Substituting the remaining field components into the 
linearized Hamiltonians yields equality up to a sign, 

= -Vl/t ^ ^AB ^ ^KK'^j 

which shows that the Zeeman effect and the ISO coupling 
act equivalently on the lowest Landau level. We remark that 
this reasoning is only valid at this specific Landau level: For 
higher Landau levels, all the eight fields are present, and the 
two terms become inequivalent. 

2. The spin-imbalanced quantum Hall phases 



1. The Zeeman-induced weak QSH phase 

As we already mentioned, the spectrum of the spin- 
degenerate model shows Landau levels for low values of the 



flux (cf. Sec. IV A I. The Hall conductivities in the gaps are 
a-a/i^e^ /h) — . . . , —6, —2, 2, 6, . . .. If we concentrate on the 
bulk gaps around zero energy, we observe one and then three 
doubly degenerate edge states (see Fig.|5ja)). When the Zee- 
man splitting is nonzero, the two copies of these states shift 



So far, most of the studies on the QSH effect have gener- 
ally concentrated on the behavior of the system at zero energy. 
However, very interesting features also emerge at nonzero en- 
ergies in the presence of external fields. For instance, in the 
low-flux regime, a Zeeman gap is formed at the n ~ 1 Lan- 
dau level (at E/t « 0.42), as shown in Fig. |5jb). Similarly 
to the zero-energy gap, this gap also shows a difference be- 
tween the number of edge states with spin up and spin down 
components, A'^ = 1 and — 3, which may be understood 
from the values above and below the corresponding Landau 
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FIG. 5. (Color online) Comparison of dispersions at <j!> = 1/61 for (a) the spin-degenerate case, (b) the system with Zeeman effect (Az jt = 0.5) 
and (c) with ISO {tijt — 0.01). The coupling constants Az and ti have been chosen such that the (weak-QSH) gap sizes at zero energy become 
equal. 




FIG. 6. (Color online) High-energy spectra for the model without 
spin-orbit couplings (ta = ti = 0) at = 1/61. (a) spectrum for 
the spin-degenerate model (Az = 0), with quantum Hall gaps cor- 
responding to (jh = 2, 4, 6, . . . (in units of the conductivity quan- 
tum), (b) Spectrum with Zeeman effect (Az ~ 0.2t), showing the 
spin-filtered quantum Hall state with A^^ — 1 and Ni — 0. The 
other gaps show (alternatingly) quantum Hall phases (A^^ — Ni) 
and spin-imbalanced quantum Hall phases (\Nf\ = \Ni\ + 1). 



spective conductance quanta. The spin-filtered quantum Hall 
phase appears at high energies and low flux in the presence of 
Zeeman coupling, as we now explain. In the absence of Zee- 
man coupling, we observe equally-spaced Landau levels, each 
of which has Chern number —2, see Fig.|6|a). If the Zeeman 
term is present, the gap that forms between the two copies 
of the original highest-energy Landau level exhibits a spin- 
filtered quantum Hall phase, characterized by the presence of 
edge states of only one spin component; in this case A'^ = 1 
and — 0, see Fig. [6|b). The other gaps visible in Fig. |6|b) 
are spin-imbalanced and ordinary spin-degenerate quantum 
Hall gaps, alternatingly. Let us mention that the spin-filtered 
and spin-imbalanced quantum Hall phases are ubiquitous in 
systems with Zeeman-split Landau levels, e.g., the quantum 
Hall plateaus corresponding to odd-integer filling factors in 
GaAs/AlGaAs heterostructures,'^ or in HgTe quantum wells, 
which have strong Zeeman effect and are thus ideal candidates 
for observation of these phases.ES 



level of the spinless model. Again, the spin Hall conductiv- 
ity is nonzero, = — 2e/47r, but the difference with the 
weak quantum spin Hall gap is that the charge Hall conduc- 
tivity is also nonzero, cth = 4e^//i, as can be deduced from 
Fig. ^h). Furthermore, the edge states in this gap all propa- 
gate in the same direction (i.e., they are chiral), thus providing 
robustness against disorder Here, we refer to this phase as the 
spin-imbalanced quantum Hall phase!^ 

A special instance of the spin-imbalanced quantum Hall 
phase may occur if edge states of one of the spin-components 
are absent, i.e., either = or A^^ = 0, in which case 
one speaks about the spin-filtered quantum Hall phase. In this 
state, the magnitudes of the charge Hall and spin Hall cur- 
rents are therefore equal when expressed in units of their re- 



D. Rashba spin-orbit coupling 

The RSO coupling differs in an essential way from the ISO 
coupling and the Zeeman effect, in the sense that the hopping 
involves a spin flip. In other words, the spin matrices involved 
in the Rashba hopping are off-diagonal (i.e., a.^. and cjy), 
whereas the other terms in the Hamiltonian involve diagonal 
spin matrices. This has a profound effect on the spin structure 
of the system: If the RSO interaction is nonzero, mixing oc- 
curs between the spin components and we find that the spin 
orientation of the eigenstates is generally site-dependent. The 
spin direction is confined to the two-dimensional plane and is 
perpendicular to the momentum: Aside from the out-of-plane 
component (z-direction), the component perpendicular to the 
cylinder edge may also be nonzero, while the component par- 
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model discussed in Sec. Ill A will cause a spin-splitting of 
some bulk bands in the Hofstadter butterfly, see Fig. |7] In 
particular, at low energies, the Rashba term splits the fourfold- 
degenerate Landau levels into twofold-degenerate Landau lev- 
els. In the limit — > 0, we also note that the spectrum indeed 
remains gapless and particle-hole symmetric. However, this is 
not true for other Dirac regimes in the spectrum. For instance, 
at = 7/15 and E/t ss 1.7 (see the inset of Fig.|7]), the RSO 
coupling splits the central Landau level, and the resulting gap 
exhibits two counterpropagating edge modes on each edge. 
However, the spin Hall conductivity is undefined because the 
spin direction varies as a function of the Fermi energy. Never- 
theless, at each energy inside the gap the spins of the two edge 
states are approximately equal, so that the spin current nearly 
vanishes. Besides, the charge Hall conductivity is exactly zero 
when the Fermi energy is located in this gap. 



FIG. 7. (Color online) Hofstadter butterfly spectrum (with RSO cou- 
pling, t-R/t — 0.2). The inset shows a magnification of the region 
surrounded by the dashed lines. In the inset, the dashed line indi- 
cates the flux value (ji — 7/15 discussed in the text. The colors and 
numbers indicate the quantized Hall conductivity in units of /h. 



allel to the edge always vanishes. In other words, the helicity 
{o'x, o'yjO'z) ■ fik of the (edge) states remains zero. This effect 
is caused by the nature of the Rashba coupling in the Hamil- 
tonian ([T]), which involves the outer product of the spin vector 
(T with the bond vectors djk, so that the spin is always per- 
pendicular to the hopping direction. Because the edge-state 
spins are no longer in the up or down state, the Chern numbers 
for the respective spin components are no longer well-defined. 
Thus, the spin conductivity no longer takes quantized values. 
Nevertheless, the spin Chern number remains a well-defined 
integer-valued topological invariant,**' which can be used as a 
tool to distinguish between trivial and nontrivial gaps. 

The second difference with the ISO coupling is that the 
zero-flux spectrum remains gapless and the spin degener- 
acy is lifted by the RSO coupling (considering tR ^ and 
ti ^ Xz = 0). Indeed, upon setting tj^ ^ 0, the two Dirac 
cones from the spinless model (at the special points K and K') 
are broken into four cones each: A central isotropic cone at K 
or K' , and three anisotropic "satellit e" co nes around it. This 
effect is known as trigonal warping-I^^ISil Between the main 
cones and the satellite ones, there is a Van Hove singularity 
at low energy, which scales as ±(tR/t)^ for small t-R jt. The 
Berry phases associated to the main and satellite cones are 
— TT and TT, respectively, which add up to the Berry phase 27r 
of the Dirac cones in the spin-degenerate model (tt for each 
spin component). For energies |£'/t| < (in/i)"^, the low- 
flux Landau level spectrum is characterised by two different 
sets of Landau levels with Chern numbers 2 and 6. Outside 
this regime, there are only twofold-degenerate Landau levels, 
because the satellite cones can no longer be resolved. Around 
\Elt\ w (tri /t)^ there is a crossover regime where the sixfold- 
degenerate bands split into three, each one being twofold de- 
generate. 

Adding a nonzero Rashba term to the spin-degenerate 



The RSO coupling has drastic effects on the edge states, 
as their spins no longer align perpendicularly to the sample 
plane, but get an additional in-plane component perpendicular 
to the edge. Besides, the direction of the edge-state spins de- 
pends on the Fermi energy. This observation distinguishes the 
Rashba effect from the Zeeman effect in a tilted field, which 
puts all edge-state spins in the same direction. In addition, the 
in-plane component is opposite for edge states at the oppo- 
site edges, whereas the perpendicular component is the same. 
Importantly, due to the dependence of the spin direction on 
the Fermi energy, the latter can be conveniently used for spin 
manipulations (in addition to the amplitude of the Rashba cou- 
pling itself). 



We illustrate this phenomenon around the Rashba-split bulk 
band at = 7/15 and Ejt « 1.7. In Fig.jsfa), we show that 
in the limit of ^r 0, the eigenstates are in-plane (indicated 
as "left" and "right" in the figure, meaning that the spin di- 
rection is ±y). This result shows that for infinitesimal RSO 
coupling, it is more natural to decompose the states in terms 
of the eigenstates of the Pauli matrix Gy than those of a^. If 
the Rashba coupling is then increased, the gap opens, and in 
the vicinity of the bulk bands, the edge states tend to rotate to 
the vertical direction, see Fig. |8jb). In addition, we observe 
that not all spins have the same length. This phenomenon oc- 
curs because the spins displayed in this plot represent the ex- 
pectation values of the spin components (see the Appendix). 
Since the spin direction depends on the lattice position (y co- 
ordinate), the length of this expectation value may be less than 
unity. 



We expect that the canting of the spins is most easily ob- 
served in gaps with a single edge state on each edge. For 
instance, in Ref. 9, it has been demonstrated that the RSO 
coupling cants the spin in a spin-filtered gap generated by the 
Zeeman effect. In this setting, the spin textures may be con- 
trolled at will by tuning the coupling parameters and the Fermi 
energy. 
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FIG. 8. (Color online) Dispersions at = 7/15, showing the rotation of spin directions as the Rashba coupling tR is increased. 



Comparison between the low-flux limit < 
"half-flux" regime (j)'^l/2 



and the 



Az ^ 



\z/t = 0.01 




FIG. 9. (Color online) Spin splitting in an example at flux — 7/15 
and near E/t = 1.7t. Zeeman effect: (a) Az — >■ and (b) Az = 
O.Olt; ISO coupling: (c) ti and (d) tj = O.Olt; RSO coupling: 
(e) tR and (f) tn = 0.2t. 



As shown before in the low-flux limit </> — > 0, the ISO 
coupling and the Zeeman effect affect the zero-energy modes 
equivalently in the sense that they both open a weak QSH 
gap. Indeed, it has been shown that for this mode, the ISO 
coupling and Zeeman effect are formally equivalent."^"* On the 
other hand, the RSO coupling does not open a gap at zero en- 
ergy. This fact hampers the comparison of the RSO coupling 
to the aforementioned two terms in this energy regime. 

From the Hofstadter butterfly shown in Fig.[T] one observes 
that the Landau-level structure in the low-flux, low-energy 
limit is not unique. In fact, similar structures appear wherever 
the shape of the dispersion is characterized by Dirac cones. 
For example, at = 1/2 and for energies close to E/t = 
±-\/3, the dispersion can be approximated by a Dirac disper- 
sion of the form \E/t\ = y3±y||k|- i\/3|k|2 + 0(|k|3), 

where k is the momentum relative to the position of one of the 
two Dirac cones. The linear term is responsible for the square- 
root behavior of the (fourfold-degenerate) Landau levels for 
flux values close to <j) ~ 1/2. The quadratic term causes the 
asymmetry of the Landau level spectrum with respect to the 
energy \E/t\ = \/3. Furthermore, we observe that the energy 
\EhL.o/t\ ~ ^3 - (7r/3) 10- 1/2 1 of the central Landau level 
is linear in the flux rather than constant. 

The approximate linear dependence of the spectrum i?(k) 
at<f> = 1/2 and \E/t\ w -\/3 suggests that the Zeeman and ISO 
terms should behave similarly in this regime, in direct analogy 
to the case (f>,E ^ discussed in Sec. IV C 1 Although the 
energy of the lowest Landau level is no longer constant in cj), 
the approximation holds, thus leading to a similar behavior 
with regard to the opening of the gaps, as shown in Figs.[9|a- 
d) for the flux value </> = 7/15 which lies close to (p = 1/2. 
The Zeeman effect and ISO coupling split the central Landau 
level at E/t w 1.7 and open up a weak QSH gap, similar to the 
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zero mode in the low-flux regime. As discussed in Sec. IV D 
the RSO coupling splits the central Landau level, unlike the 
case for — >^ 0, see Figs. |9|e,f). Although there are edge 
states inside the gap, it is a trivial state, because the spin Hall 
conductivity approximately vanishes. This comparison shows 
an important difference between Zeeman coupling and ISO 
coupling on one hand, and the Rashba coupling on the other 
hand: The former two create the weak QSH phase, while the 
latter does not. In the following section, we show that in the 
presence of Zeeman or ISO coupling, the Rashba coupling 
tends to destroy the weak QSH phases generated by the former 
effects. 



V. COMPETITION AND PHASE TRANSITIONS 

After analyzing the effects of the three terms Hi, Hj^, 
and Hz in Hamiltonian ([T]l independently, a natural question 
arises: what effects emerge when these terms are combined? 
Here, we study the competition featuring two terms, by tun- 
ing the ratio between their corresponding coupling constants. 
Tuning these coupling constants generally leads to opening 
and closing of gaps, so that the topological nature of the gaps 
may change. In that case, we deal with topological phase tran- 
sitions. We recall that Kane and Mele"^ have studied the phase 
transition from a (time-reversal symmetric) QSH state to a 
trivial state, by increasing the ratio ^R/ti. Here, we extend 
this study and we explore the more exotic phase transitions 
that are realized in the presence of a magnetic field and the 
Zeeman effect. 



A. Phase transitions driven by the ISO coupling 

Let us analyze the interplay between the ISO coupling and 
the magnetic field by investigating a gap at a fixed flux value, 
in absence of the Zeeman effect and of the Rashba coupling. 
By tuning the value of the ISO coupling amplitude ti, we 
can close the gaps and obtain different topological phases 
at both sides of the transition. In Fig. [TO] we have shown 
the dispersions for = 1/3, for values of ti below, at, and 
above the values where the gap around E/t — 2.3 closes. 
Around this energy, there are three bulk bands: Two (over- 
lapping) ones above and one below the gap. For values of 
ti < ti,o ~ 0.4285t, below the critical value, we observe the 
weak QSH phase, i.e., two counterpropagating edge modes 
with opposite spins on each edge, = 1 and = — 1. The 
spin-down pair of edge modes connects the top and bottom 
bulk bands visible in this plot, and the spin-up pair connects 
the middle band to a lower band. For ti — o, the gap dis- 
appears: The bulk bands touch each other in three points at 
E/t ss 2.3, and two of these points are connected by the spin- 
down edge states in the middle of the figure. If the value of tj 
is increased further, a gap opens again, and a new pair of edge 
states (with spin down) appears. Moreover, the existing spin- 
down edge states have inverted their direction of propagation. 
Thus, for ti > ti 0, we have a spin-imbalanced quantum Hall 
state with A^i^ = 1 and iV^ = 2. Note that only the spin-down 



states are modified in this phase transition, while the spin-up 
states remain the same. 

In Ref.[9] this phase transition was investigated in the pres- 
ence of a nonzero Zeeman coupling in addition to the ISO 
coupling. Here, we emphasize that the Zeeman effect is ac- 
tually superfluous and that this transition is driven exclusively 
by the ISO coupling. Variation of Az has no qualitative ef- 
fect on the transition, but only shifts the energies of the bands 
and gaps. Within a large range of values for Az, we obtain an 
identical phase transition as in Fig. 10 The critical value o 
of the ISO strength is independent ofAz. 

The phase transition presented here differs in a fundamen- 
tal way from a crossing of Landau levels. In the latter case, 
the Chem numbers associated to the Landau levels do not 
change at the transition, and as a consequence the number of 
edge states inside the gap between them remains unmodified 
as well. For the phase transition presented here, the number 
of edge states and consequently the Chern numbers do change 
at the transition. Focusing on the spin-down edge states, we 



observe that N 



-1 for ti < ti^o and iV; = 2 for > <i 



a difference of AA'^ — 3. The Chern numbers of the bands 
below and above the gap are 1 and 1 in the former case and 
—2 and 4 in the latter. The sum of these Chern numbers is un- 
changed, as required by the bulk-boundary correspondence, 
Eq.(|5j. 

The question arises as to whether the difference in the num- 
ber of edge states at both sides of the transition may be pre- 
dicted. A priori, this difference cannot be predicted since it is 
the result of a complicated interplay of the ISO coupling and 
the magnetic flux. However, the difference is always a multi- 
ple of the denominator q of the flux (j> = p/q. This result may 
be understood from the fact that at a given flux value (j) = p/q, 
the Chem numbers of all the Hofstadter bands obey C„ = c 
(mod q). Here, c is the modular multiplicative inverse of p 
modulo q, defined as the unique integer c (0 < c < g) such 
that cp=l (mod g).— Consequently, the difference between 
two possible values of the Chern numbers is always a multiple 
of q, which proves the aforementioned claim. This property 
may also be understood from the superlattice structure: The 
magnetic field produces q copies of the original unit cell, and 
each of them reacts in the same way to the perturbation lead- 
ing to the phase transition. 

In Fig.[TT| we show the phase diagram as a function of the 
ISO coupling ti and the Fermi energy Ep with the flux value 
(/i = 1/3 held fixed. From this figure, we can easily read off 
the values of <i and the energies E-p where the gaps close and 
the topological phase transitions occur. A few points must 
be noticed: First of all, we clearly identify the weak-QSH 
phase at half-filling for ti/t = 0.4. The (spin-degenerate) QH 
phases are progressively destroyed for ti/t < 0.3. For h/t > 
0.3, we get spin-filtered QH [e.g., iN^,Ni) = (-1,0)], 
spin-imbalanced QH [e.g., (iV^,iV;) = (+l,+2)] and the 
weak QSH phase. In addition, there are also more exotic 
phases where the edge states are neither chiral nor helical, 
e.g., {N^,Ni) — (+1,-2). The spin-degenerate phases do 
not arise for ti/t > 0.3. Secondly, the gaps are large (of the 
order of t) in wide regions of the parameters, which is favor- 
able for the detection of the transitions in cold-atoms experi- 
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FIG. 10. (Color online) Spectra for = 1/3, with (a) ti/t = 0.35, (b) U/t = 0.4285, (c) ti/t = 0.5. In (a), we observe a (weak) QSH gap, 
where TV-j- = —Ni — 1. (b) By increasing the ISO coupling constant, the bulk gap closes for ti/t ~ 0.4285. (c) If ti is increased further, 
the gap reopens, but with a different number of edge states, namely A''^ — 1 and N_i = 2, so that it enters a spin-imbalanced phase. Thus, the 
system undergoes a topological phase transition from a helical phase (weak QSH) to a chiral phase (spin-imbalanced QH). 



EF/t 




FIG. II. (Color online) Phase diaj 
pling ti and the Fermi energy Ef 
In the shaded regions, the system 
gray) and blue (light gray) disting 
ponents, respectively. In the white 
the edges conduct. There, the pair 
of edge states A'^ and A''^ for spin 



;ram as a function of the ISO cou- 
for the fixed flux value (j) — 1 /3. 
is metallic. The colors red (dark 
;uish the spin-up and -down com- 
regions, the bulk is insulating, and 
s of numbers indicate the number 
up and down, respectively. 



that the Zeeman effect alone does not induce the same type 
of phase transitions as the ISO coupling. The Zeeman effect 
can only close and open gaps between bands of different spin 
components, and does therefore not modify the Chem num- 
bers of the bands, whereas the ISO coupling can also close 
gaps between two bands of the same spin component. Nev- 
ertheless, tuning the Zeeman coupling strength Az allows one 
to modify the nature of the phase transitions driven by ISO 
coupling. 

B. Intrinsic spin-orbit coupling and Zeeman effect for 4> ^ 

If the Zeeman effect and ISO coupling are present simul- 
taneously, the interplay between the two terms is governed 
by the value of the flux cj), because the amplitude of the for- 
mer is linear in while that of the latter is constant. Thus, 
in the butterfly spectrum illustrated in Fig. [12] we observe 
the ISO-dominated regime at low flux (cf Fig. |3]l, and the 
Zeeman-dominated regime at high flux (cf. Fig. [4]). The 
crossover between these two regimes takes place where the 
two terms are comparable in strength, approximately where 
intfiXz ^ GVSti, indicated by the vertical dashed line in 
Fig. 



12 



at w 0.165. We remark that this crossover effect 
is present more generally in the situation where the Zeeman 
term competes with any other term with a constant amplitude: 
Due to its linear ^-dependence, the Zeeman term will always 
dominate any other one at sufficiently high values of (p. Natu- 
rally, this phenomenon only occurs in systems where the Zee- 
man splitting depends linearly on the magnetic flux (p. In the 
presence of magnetic impurities, the Zeeman term may ex- 
hibit a nonlinear dependence on the applied magnetic flux, as 
in Fig. 1 1 up (down) by a fixed amount of energy. We note discussed in Ref. 49 for Mn-doped HgTe quantum wells, but 



ments. In addition, the property that the number of edge states 
changes by a multiple of 3 at each phase transition is clearly 
observed in Fig. 
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Finally, the energies where the phase 
transitions occur may be shifted by tuning the strength of the 
Zeeman effect. Since the spin-up and spin-down components 
are uncoupled, the Zeeman effect will shift all red (blue) areas 
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FIG. 12. (Color online) Hofstadter butterfly spectrum with ISO cou- 
pling {ti/t — 0.1) and Zeeman effect (\z/t — 0.5). the colors and 
numbers inside the bulk gaps indicate the quantized Hall conductiv- 
ity in units of /h. The vertical dashed line indicates the flux value 
where the two effects are comparable in strength. At this value, the 
crossover from the ISO-dominated (low </>) to the Zeeman-dominated 
(high (fi) regime takes place. 



this subject is beyond the scope of this paper. 

Although the Zeeman effect dominates over the ISO cou- 
pling for (p ^ 0.165, the influence of the ISO coupling is still 
visible in the high-flux spectrum: At = 1/2 and </> = 1, 
there are gaps at E/t 



12 1, which are absent in 



3 (see Fig. 

the butterfly without the ISO coupling (see Fig. |4|. We note 
that the dispersions around these points are of Dirac type for 
ti = 0, consequently, the ISO coupling opens these gaps in a 
way analogous to the gaps at zero flux and zero energy. How- 
ever, the gaps at (f> — 1/2 and (f) = 1 are spin-filtered QH 
phases and not QSH, due to the subtle competition between 
the ISO coupling and the TRS-breaking effects. 



C. Rashba spin-orbit coupling and exchange term at = 

In two recent works,'^it has been shown that the combina- 
tion of RSO and an exchange field, similar in structure to the 
Zeeman term, leads to chiral edge states at zero energy. In ad- 
dition, the problem has been investigated for a sodium/lithium 



iridate model which involves NN and NNN hopping terms.'Sl 
In the latter and more complex setup, nontrivial QH gaps ap- 
pear also away from half filling. Here, we consider the com- 
bined effects of the Zeeman and Rashba couplings, away from 
half-filling and in the absence of ISO coupling. We demon- 
strate that several QH gaps appear in the absence of NNN hop- 
ping terms, without including any gauge field (i.e., the Peierls 
phases): Namely, non-trivial QH phases are produced in our 
model by setting (p = ti ~ Q, ts. ?^ an d by applying a 



(8) 



flux-independent Zeeman, or exchange ,^1-^ term 



which differs from by the fact that the strength ^ox no 
longer depends on the flux (p (which we now set to zero). We 
stress that it is the simple association of a constant exchange 
term with a NN Rashba hopping term that leads to non-trivial 
topological phases, at half-filling but also at E-p « gcx, as 
discussed in detail below. 

When the exchange term is added to the tight-binding 
model with strength < gcx < 3i, the energy bands asso- 
ciated to the spin-up and spin-down components are shifted in 
opposite directions but overlap at half-filling. In this situation. 



which is illustrated in Fig. 13 a), the system is a semi-metal 
for Ep = ±(7cx and a metal otherwise. When the Rashba 
coupling is turned on, bulk gaps open not only at half-filling 
E = but also at E = ±.gcx (cf. Fig.[T3|b)). We have com- 
puted the edge-state structures as well as the topological in- 
dices associated to the gaps: we obtain that all these gaps are 
related to non-trivial QH phases. More precisely, we find that 
the gaps at E = ±.gcx host a single edge state per edge, the 
spin orientation of which is given by (ctz) « ±1: these gaps 
coiTespond to spin-filtered QH phases, with equal (resp. op- 
posite) charge (resp. spin) Hall conductivities. We observe 
that the spin orientation of these isolated edge states is nearly 
vertical, which is due to the fact that the exchange term is im- 
portant in this high-energy regime. At half-filling, the bulk 
gap hosts two edge states per edge, with same velocity and 
quasi-horizontal spin orientation. This bulk gap coiTesponds 
to a topological QH phase with cth = 2 and a non-trivial spin 
orientation. Therefore the spin orientation of the edge states 
highly depends on the Fermi energy. We note that the gap 
openings occur for large exchange coupling, but for arbitrar- 
ily low Rashba coupling. 

When the exchange strength > 3t, the energy bands 
associated to the spin-up and spin-down components are well 
separated: a trivial gap opens at half-filling, while the Dirac 
points are shifted to the energies E — ztg^x, cf. Fig. [T3|c). 
In this situation, the Rashba coupling still opens two non- 
trivial QH phases at Ep = ±gcx, cf. Fig. [T3|d). Again, 
these phases are characterized by spin-filtered edge states with 
nearly vertical spin orientation. Therefore, for sufficiently 
large exchange coupling, it is necessary to set the Fermi en- 
ergy away from half-filling in order to detect nontrivial topo- 
logical phases. This behavior is summarized in the phase di- 
agram of Fig. 13 e): For nonzero Rashba coupling tjx, the 
Hall conductivity at zero energy is cth — ±2 for \gcx/t\ < 3 
and zero otherwise. Thus, topological phase transitions of the 



13 



t^/t = 0, 3ex/t = 2 




tn/t = l,Pox/t = 2 




(a) 


III! 

k 

tn /t = 0, gex/t = 4 


(b) 


tuft 


' ' ' 'fe 
= l,SrexA = 4 


E/t 
5 - 




5 - 






- 




- 
















-5 - 




-5 - 












(c) 


1 1 1 1 , 
k 


(d) 


III! 

k 



(+1,0, +1) 


(+l,-2,+l) 


(-l,+2,-l) 


(-1,0,-1) 


(+1.0, +1) 

1 


(+l,-2,+l) 
1 


(-l,+2,-l) 
1 


(-1,0,-1) 
1 



Xz/t = 0.5, tn/t = 0.02 




0.05 



-0.05 



FIG. 14. (Color online) Spectram for the system with nonzero Zee- 
man (\z/t = 0.5) and RSO coupling (tB./t = 0.02), for (/> = 1/61. 
The inset shows a zoom of the region where a gap opens between 
the edge states at zero energy. For comparison with the spectrum in 
absence of Rashba coupling, refer to Fig.[5|b). 



and [T3|d). This competition between QH gaps of different 
origins leads to complex and rich quantum phase transitions. 



, , -4 -2 2,4 

(e) gex/t 



FIG. 13. (Color online) (a-d) Energy spectra E — E{k) for 
(j) = ti = 0. The Zeeman (exchange) and Rashba coupling strength 
are indicated, (e) Phase diagram in the parameter space defined by 
the coupling strengths gox and fR. The number triplets indicate the 
charge Hall conductivity in the lower, middle, and upper gaps, re- 
spectively. The solid lines indicate the phase transitions for the gap 
at zero energy. At the dashed lines, the lower and upper gaps close 
and the trivial gap at zero energy remains open [as in (c)], and at the 
dotted lines, all three gaps close [as in (a)]. At the dashed and dotted 
lines, the charge Hall conductivities in the gaps do not change. 



zero-energy gap take place at the gap closings at gex/t = ±3. 
The lower and upper gaps (at Ep w ±5cx) always exhibit a 
Hall conductivity of ±1, where the sign depends on that of 
gcx/t but not on that of tn/t. Inverting the sign of gcx inverts 
the directions of the edge currents, as well as the z component 
of the spins. Inverting the sign of % inverts the in-plane com- 
ponent of the spins. Continuous variation of in through zero 
changes the spin-direction continuously, where the spins tend 
to the z direction in the limit t^ /t 0. Thus, the gap closing 
att^ — does not constitute a topological phase transition. 

Let us emphasize the fact that the chiral QH phases pre- 
sented here are not the consequence of the external mag- 
netic flux (i.e., Peierls phases), since they are produced when 
(j) = 0, and solely in the presence of the exchange and Rashba 
terms. When 7^ 0, new QH gaps open in the system, altering 
the exchange-Rashba-induced phases presented in Figs.[T3|b) 



D. Zeeman effect and Rashba spin orbit coupling for 0/0 



In the low-energy regime, where the (weak) QSH phase is 
generated by either ISO coupling or the Zeeman effect, the ad- 
ditional inclusion of the Rashba coupling has several profound 
effects on the physics. First of all, the RSO coupling tends to 
decrease the size of the gap (cf. Ref. 4), eventually completely 
closing it. Secondly, the combination of RSO coupling with 
ISO coupling and/or Zeeman effect destroys the particle-hole 
symmetry of the spectrum. Finally, at finite magnetic field, 
the RSO coupling opens up a gap (i.e., an avoided crossing) 
between the edge states at their crossing around zero energy. 
Around this crossing, the spin direction rotates from down to 



up (or vice versa), as illustrated in Fig. 14 The size of this gap 
is approximately linear in (p and in <r (in the limit — > and 
tft — > 0). Thus, the bulk gap is trivial in the sense that no edge 
states cross it from one bulk band to the other. However, in the 
limit where the edge-state gap is small, one would still observe 
a weak QSH state. At finite temperatures, the thermal energy 
can "bridge" the edge-state gap, so that the charge carriers 
may tunnel from one edge state to another As the edge-state 
gap size is increased, the backscattering is enhanced, due to 
the decrease in the amount of tunneling across the edge-state 
gap. Here, we expect that the value of the spin Hall conductiv- 
ity will deviate from its quantized value ±2e/47r. The amount 
of deviation increases when the gap size increases. 
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VI. EXPERIMENTAL REALIZATIONS 



that causes reentrant behavior of the topological phases.' 



ED 



In this section, we will discuss the several experimental 
approaches that allow for the observation of the topological 
phases and phase transitions discussed in Sees. IV and[V] For 
each of them, we present the opportunities and challenges, and 
we discuss the range of parameters that can be probed. 



A. Condensed-matter systems 

Condensed-matter systems typically have a lattice spac- 
ing in the order of 1 A, so that one flux quantum per unit 
cell {(p = 1) corresponds to a magnetic field in the order of 
10^-10^ T. Thus, only the low-flux limit is relevant with re- 
spect to realistic condensed-matter realizations of our model. 
The prime example of a two-dimensional system with a hon- 
eycomb lattice is graphene, which has received rnuch atten- 
tion both at the theoretical and experimental levels.'^^'^In the 
presence of a magnetic field, graphene shows an anomalous 
quantum Hall effect, where the Landau-level energies'*^ de- 
pend on the magnetic field as y/<l>, showing a radically dif- 
ferent behavior than in semiconductor heterostructures (e.g., 
GaAs/AlGaAs). 

Besides, the effects of spin-orbit coupling have been ex- 
tensively studied for graphene, leading to the concept of Z2 
topological insulators.'* However, recent analysis has shown 
that graphene exhibits only a very weak I SO coupling (ti/t ^ 
10^^-10"^, where t - 2.8 meV),EISl which unfortunately 
prevents the observation of the QSH state in this mate- 
rial. Nevertheless, the ISO interactions may be enhanced by 
putting heavy adatoms on top of the graphene surface,E2which 
would open up a route to build a topological heterojunction.D^ 
The Zeeman effect is of appreciable order of magnitude, g ^ 
2, which corresponds to gaps of 100 K at a magnetic field of 
30 T.l^ Values of the Rashba coupling as high as <r ^ O.lt 
have been reported for graphene on a Ni(lll) substrate. Ad- 
sorption of metal atoms on the graphene surface increases this 
coupling as well.^'^ Alternatively, high values of the Rashba 
coupling (accompanied by the spin re-orientation at the edges) 
may be induced by curving graphene nanoribbons.'^ Besides, 
it has been suggested that large Rashba coupling leads to topo- 
logical phases in bilayer graphene.''' In addition, large effec- 
tive pseudomagnetic fields, greater than 300 T, have been re- 
cently realized in strain-induced graphene nanobubbles.l^ In 
summary, the observation of the topological phases presented 
here remains elusive in graphene, but should not be defini- 
tively ruled out. 

The QSH effect has nevertheless been observed in other 
two-dimensional compounds that have considerably larger 
values of ij, such as Hg(Cd)Te quantum wells,!^ w here values 
of the order ti ^ 1-10 meV have been reported .'^^'^ Addi- 
tionally, this material exhibits a strong Zeeman effect as well 
(g ~ 20).^ The spin-filtered QH states have already been 
observed in these systems at low fields. Hg(Cd)Te quan- 
tum wells doped with magnetic ions (e.g., Mn) exhibit an ex- 
change interaction leading to the quantum anomalous Hall ef- 
fect at zero magnetic field, and to a nonlinear Zeeman effect 



B. Artificial honeycomb lattices 

The experimental limitation of the low flux values can be 
overcome by studying systems with a larger lattice spacing. 
Recently, several approaches have led to the engineering of ar- 
tificial honeycomb lattices. First of all, arrays of quantum dots 
on a GaAs, arranged in a honeycomb pattern, have been engi- 
neered to simulate actual graphene.'^ ^'^ The lattice spacing in 
these systems is of the order of 100 nm, so that high flux val- 
ues {4> ^ 1) are obtained already with small magnetic fields. 
However, this large lattice constant leads to a small hopping 
amplitude, which has hampered the observation of the Hofs- 
tadter butterfly until now. Patterned superlattices able to probe 
the Hofstadter regime with usual magnetic field intensities 
have been realized in the past, but for square geometries.t^^^ 
To probe the physics described here, one needs a superlattice 
with a smaller lattice constant, grown on a substrate that could 
lead to the generation of a strong Rashba interaction, for ex- 
ample. In this case, chiral topological states could be probed 



also away from zero energy, as we discuss in Sec. IV E 



Secondly, artificial graphene lattices have been created by 
manipulating CO molecules with the tip of a scanning tun- 
neling microscope and arranging them in a triangular lattice 
on a Cu(l 1 1) substrate.'"* These molecules repel the electrons 
on the Cu surface, thus creating a honeycomb array of sites 
between which the electrons hop. Measurements of the den- 
sity of states show that the electrons indeed exhibit a Dirac 
dispersion. The typical lattice constants of these systems are 
approximately 7 times larger than in real graphene, thus al- 
lowing for flux values ^ 50 times larger than in graphene, 
with reasonable magnetic fields. In addition, effective mag- 
netic fields up to 60 T can be induced by triaxial strain of the 
artificial lattice. It can be speculated that by using a heavier 
metal as substrate, it may be possible to engineer such a lattice 
with strong spin-orbit coupling, leading to the QSH phase. 



C. Ultracold atoms in optical lattices 

Although the parameters of the Hamiltonian ([T]i can be con- 
trolled to some extent in condensed-matter systems, the pos- 
sibility to engineer the physics discussed here with ultracold 
atoms in optical lattices is very appealing. The advantage of 
the ultracold-atom realization of this model is the flexibility to 
control the parameters separately across a large range, com- 
pared to condensed-matter systems where these parameters 
are gener ally fi xed by the material properties and the sample 
geometry.l^SE] A second major advantage of ultracold atoms 
over condensed-matter systems is that the magnetic fields are 
produced synthetically, i.e., the Peierls phases are not due to 
a real magnetic field, but induced by external fields acting on 
atomic internal states. In this way, large effective magnetic 
flux can be produced (i.e. (p ^ 1), which allows the system 
to be studied beyond the low-flux limit. ^' The spin-orbit cou- 
plings can also be synthesized in optical lattices using simi- 
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lar methods!^ Moreover, in these configurations, the Zeeman 
coupling can be controlled independentl y, as i t is not related to 
the synthetic magnetic fieldJ^ (see Sec. V C where we study 



the extreme situation where the effect of the Zeeman coupling 
is investigated in the absence of synthetic magnetic field, in 
which case this term is called an "exchange term"). Another 
important feature offered by optical-lattice setups is the pos- 
sibility to control the interactions between the particles: by 
exploiting Feshbach resonances,'23 it is possible to reach the 
non-interacting regime and from there on, to include interac- 
tions in a controlled way to explore the robustness of topolog- 
ical order against interactions. Let us mention that these sys- 
tems are usually free of disorder, although disorder can also 
be generated and tailored at will in optical lattices.l^l These 
outstanding features allowcus to study the TRS-broken QSH 
phase, namely the survival of the QSH phase in the pres- 
ence of a TRS-breaking perturbation .121 Indeed, contrarily to 
condensed-matter experiments, the absence of (magnetic) dis- 
order can stabilize the helical (counterpropagating) edge states 
even when TRS is broken. Let us also comment on the fact 
that all the topological phases and edge-state excitations dis- 
cussed here rely on the existence of robust bulk energy gaps. 
Therefore, when the temperature of the system approaches the 
size of the topological gaps of interest, one expects the effects 
to be strongly affected, and eventually to disappear For the 
ultracold temperatures achieved in cold-atom laboratories, the 
gaps of interest should typically be of the order A ^ ^, where t 
is the tunneling amplitude. Therefore, we anticipate that only 
the largest topological gaps presented in our figures should be 
considered from the cold-atom point of view. 

The quantum phase transitions presented in this work could 
be equally characterized either by a change in the topolog- 
ical invariants associated to the bulk gaps of interest, or 
by the modification of the edge-state structures. Although 
cold atoms do not offer the possibility to directly measure 
transport coefficients, several methods have been proposed 
to detect the topological invariants related to the quantized 
conductivities.'^''^'* These proposals are based on the fact that 
non-trivial band structures modify the atomic densities, lead- 
ing to several signatures that can be directly seen in time-of- 
flight or density measurements. On the other hand, it is possi- 
ble to directly probe edge state structures using Bragg spec- 
troscopy techniques. '^^"^^ In a recent work,'^ it was shown 
that a state-dependent light probe, focused on the edge of 
the cloud and transferring angular momentum to the atoms, 
clearly identifies the presence of chiral edge states in a QH op- 
tical lattice. This method allows us to project the topological 
edge states on a dark background, making them visible in in 
situ density measurements. In principle, this efficient scheme 
could be generalized to spin- 1/2 systems in order to detect and 
distinguish between chiral and helical edge structures. 



VII. CONCLUSIONS 

We have shown an in-depth analysis of the edge-state struc- 
tures generated by the combination of the ISO coupling, the 
RSO coupling, and the Zeeman effect, together with a uniform 



and perpendicular magnetic field. In Sec. IV we have studied 
each term individually in order to pinpoint their effect on the 
quantum Hall states generated by the magnetic field. In the 
low-field limit, the ISO coupling and Zeeman effect are for- 
mally equivalent, and give rise to the weak QSH phase. Away 
from zero energy, we observe spin-filtered, spin-imbalanced, 
and ordinary quantum Hall phases. The RSO coupling is 
different in the sense that it involves spin-flip terms, which 
causes the spin states to acquire an in-plane component. In ad- 
dition, the spin direction inside the gap depends on the Fermi 
energy, and consequently the spin Hall conductivity of the sys- 
tem is ill-defined. In the case of high magnetic fields, we ob- 
serve Dirac-type regimes away from zero energy. Here, the 
Zeeman and ISO coupling open weak QSH gaps in a way 
analogous to the low-flux limit, whereas the RSO coupling 
opens a trivial gap (unlike the case at low flux, where the 
spectrum remains gapless). In addition, the topological phase 
transitions reported in an earlier worl^have been analyzed in 
more detail here. 

A combination of the spin-orbit terms and the Zeeman ef- 
fect leads to subtle competitions, several examples of which 
have been discussed in Sec.|V] We have shown that, although 
the Zeeman effect and ISO coupling both give rise to the QSH 
phase, they do not always reinforce each other, which has been 
illustrated by the size of the parameter regime of the weak 
QSH phase in the low-flux limit. At high flux, the Zeeman 
effect dominates over the ISO coupling, but the effects of the 
latter are still well visible in the spectrum. Interestingly, a 
combination of the Rashba term with an exchange coupling 
leads to nontrivial topological phases at zero magnetic field, 
away from half-filling. Finally, inclusion of the Rashba cou- 
pling in the case where we have a weak QSH gap generated 
by the Zeeman effect opens a gap at zero energy between the 
edge states, which provides an illustration of the breakdown of 
the QSH effect in the absence of TRS and spin conservation. 

This work emphasizes the various topological phase tran- 
sitions and edge-state structures that could be observed in a 
wide range of materials and quantum emulators. Although 
we focused our analysis on the honeycomb lattice, where the 
spin-orbit couplings and the magnetic field configurations are 
expressed by specific tunneling terms and Peierls phases, we 
note that our results do not rely on this specific geometry. The 
topological phases and the phase transitions discussed here 
are universal, in the sense that they only rely on the inclu- 
sion of general ingredients, which can be found or engineered 
in many different contexts and geometries. Indeed, this uni- 
versal property can be understood in the low-energy limit, in 
which various (a priori, very different) tight-binding models 
actually reduce to the same (Dirac-type) equations. In general, 
we observe that the main ingredients needed to generate such 
physics are nontrivial gauge fields, which are coupled to the 
particles of interest (e.g., electrons in materials and atoms in 
optical lattices). In this work, we have aimed to emphasize the 
effects produced by each of these ingredients independently: 
We have investigated the effects of spin splitting, spin mixing, 
and gauge potentials leading to magnetic or Haldane-type lo- 
cal flux. We hope that our detailed analysis will deepen the 
understanding of topological insulating phases and that it will 
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FIG. 15. (Color online) (a) The honeycomb lattice. We indicated 
the NN vectors (j — 1,2,3) and the lattice vectors ai and a2. 
The sites of the sublattices A and B are indicated with solid and open 
circles, respectively. The dashed lines define boundaries of the unit 
cells, labeled as (n, m). (b) Hopping phases 0jk/2TT for the unit cell 
(n, m). The phases for the NN hopping are in green. The phases for 
the NNN hopping are in red and blue for the sublattices A and B, 
respectively. 

motivate further developments in this exciting and rapidly de- 
veloping field. 
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Appendix: Geometry and gauge 

Geometrically, the honeycomb lattice is equivalent to a dis- 
torted square lattice where each unit cell has the shape of a 
rhombus. Since there are two inequivalent sets of sites, there 
are two sites per unit cell, and hence two sublattices, labeled 
A and B. The sites, labeled as j and k in the Hamiltonian Q, 
may be uniquely identified by two integer labels n and m for 
the unit cell and A or B for the sublattice: In the following, we 
shaU write j — {n, m, t) (where n,m ^ Z and r G {A, B}) 
to uniquely identify the sites. 



The coordinates of the sites are given as follows. We choose 
the bond vectors ei, 62 and 63 as displayed in Fig. [TSja), 
where their length a = 1 is chosen as the unit of length 
throughout this text. Subsequently, the lattice vectors are 
given by ai = ei — 62 and a2 — ei — 63, so that the co- 
ordinates of the sites are 

TCn,m,A = nsLi + ma.2, rn.m,B = nsLi + ma2 + ei. 

Furthermore, the reciprocal space is spanned by the two vec- 
tors bi and b2, defined by • ~ 2TT5ij. We note that the 
results derived in this paper do not depend on the values of 
the vectors e,; the only effect of a change of the is that the 
lattice vectors and the reciprocal lattice vectors will be differ- 
ent. A convenient choice for the bond vectors is ei = (0, 0), 
62 = (—1, 0) and 63 = (0, —1), effectively putting both sites 
of the unit cell {n, m) on the coordinate (n, m). 

The phases due to the magnetic field, associated with the 
hopping, are given in terms of the gauge potential A as Oj^. — 
(e/h) J^^ A • dl, where the integral is over the line between 
sites k and j. We may employ the gauge freedom for the 
gauge potential to choose it in such a way that the phase 9jk 
for two sites j = {n,m,T) and k = {n',m',T') (NNs or 
NNNs) does not depend on the unit cell index m. The hopping 
phases resulting from such a choice are given in Fig.[T5|b) for 
hopping from the sites of the unit cell {n, m). The reader may 
check that the phase picked up from hopping around a loop 
is proportional to the area of the loop. In particular, when 
one considers one hexagon as the loop, the phase picked up is 
2Tr(f), which corresponds to an enclosed flux of cf) (in units of 
the flux quantum). Here, we note that the area of one hexagon 
is 3y/3/2, so that (j) is related to the magnetic field strength B 
as = {3V3/2)eB/h. 

In addition to the gauge choice, let us furthermore assume 
that the flux value is a rational number, (p = p/q, where p 
and q are coprime integers. Then, the phase assigned to each 
hopping is periodic in the index n with periodicity q. In other 
words, the unit cell of the superlattice (i.e., the lattice together 
with the hopping phases) is g x 1 unit cells of the original lat- 
tice. With this observation, we may invoke Bloch's theorem, 
which states that we could write the electron wave functions 

, . ^ (A.l) 

With this ansatz, the Schrodinger equation which the fields 
have to satisfy reduces to a 4(7-component matrix equa- 
tion: The degrees of freedom are all fields 'I'n — 
(V'nAt: ^riAi, V'nBt; V'«b;) in the Unit Cell of the superlat- 
tice, so that n = 1, . . . ,q, and there are two sublattice and 
two spin components per unit cell of the original lattice. This 
matrix equation is then merely an eigenvalue equation of a 
4(7 X 4q matrix, known as the Harper equation,E2EIl 
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Here, the 4 x 4 matrices I?„ = 2?„ 
cells with the same index n and from unit cells with index n + 1 to n. (Note that the Zeeman term is diagonal.) Thus, Eq. \A.2\ 
can be written equivalently as Eq. 

The 4x4 matrices 2?„ and TZn may be written in terms of the sum of four contributions, that correspond to each of the terms 
in Hamiltonian Q. The contributions for the ordinary NN hopping are given by 
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where fi = 62 — e^, £2 = 63 — 61, and £3 = ei — 62 de- 
note the NNN vectors. Here, the phase factors due to the flux 
correspond to those in Fig. [TSjb). Due to the fact that the 
phases involve fractions of 4> which are all multiples of |, the 
(/i-periodicity of the butterfly spectra is 6, rather than the peri- 
odicity of 1 which is observed in absence of the ISO coupling. 
Finally, for the RSO coupling, the hopping matrices are given 
by 
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The solution of the Harper equation yields the dispersions 
of the bulk states, which are periodic over the so-called mag- 
netic Brillouin zone, q times smaller than the original Bril- 
louin zone. Effectively, the geometry of the system is a torus, 
in this case. However, to study edge states, a cylindrical ge- 
ometry (as displayed in Fig.|2]l is more convenient. The edge- 
state spectrum is obtained from the Harper equation similar 



to Eq. ( |A.2| l, but with a matrix of dimension Aw x Aw, where 
w is the width of the ribbon in unit cells, and with the blocks 
in the top-right and bottom-left corner replaced by zeros (i.e. 
TZw = 0). In this case, the dispersion depends on just one 
component of k, and its periodicity is the reciprocal lattice 
vector 

To determine states as bulk states or edge states, we derive 
the density profile as a function of the unit cell index 

n (or equivalently, the coordinate component y). For conve- 
nience, we choose 71 = to be the center of the ribbon, so that 
n = —w, . . . ,111, where w = 2iy + 1 is the width of the ribbon. 
The expectation value of this position, (n) = J2n "-l^nP is 
used to distinguish edge states and bulk states. Edge states are 
characterized by a density profile that is sharply peaked at one 
edge, so that (n) w ±w. On the other hand, if the expectation 
value (n) is close to zero, then it means that there is a sig- 
nificant contribution of the density away from the edge; such 
a state is identified as a bulk state. In the edge state spectra 
of this paper, we have used this expectation value to color the 
states: Red and blue colors indicate the two opposite edges, 
while gray is used for the bulk states. In the same way, we use 
the spin expectation values {(Ti) — J^n '^n'^i'^n (« — x, y, z) 
to gain information about the spin states. We note that this 
expectation is only reliable in the case that the spin is constant 
(or almost constant) where the density of the state is concen- 
trated. This may not be the case in the presence of Rashba 
coupling. In the case when the spin strongly depends on posi- 
tion, the length of the vector (cr) is less than unity. 
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